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We have only rules of thumb with which to predict how a material will crystallize, chief among 
which is Ostwald's rule of stages. It states that the first phase to appear upon transformation 
of a parent phase is the one closest to it in free energy. Although sometimes upheld, the rule is 
without theoretical foundation and is not universally obeyed, highlighting the need for microscopic 
understanding of crystallization controls. Here we study in detail the crystallization pathways of a 
prototypical model of patchy particles. The range of crystallization pathways it exhibits is richer 
than can be predicted by Ostwald's rule, but a combination of simulation and analytic theory reveals 
clearly how these pathways are selected by microscopic parameters. Our results suggest strategies 
for controlling self-assembly pathways in simulation and experiment. 



I. INTRODUCTION 

Crystallization frequently happens in a 'multi-stage' 
manner, with a parent phase (e.g. a solution) first trans- 
forming into an intermediate phase (e.g. a dense liquid) 
before the stable solid emerges P-E,- One of the few 
guidelines we have for predicting when crystallization in- 
termediates will appear is Ostwald's rule of stages, which 
states that the parent phase will first transform into the 
phase closest to it in free energy |4] |5] . It is widely up- 
held. For instance, sulfur crystallizes from solution by 
first forming a dense liquid [4 . Melts [6 and aerosols [7 
also display expected precursors of the stable crystal. On 
the computer, a microscopic analog of the rule is seen: 
the freezing of polar fluids, model proteins [8 and molec- 
ular nitrogen |9j can all take place via nuclei whose com- 
position reflects that of an intermediate phase. But the 
rule has no theoretical foundation [TO] and is not univer- 
sally obeyed. Amino acid crystallization [TT], the sim- 
ulated freezing of molecular CO2 [12 and Potts model 
phase transformations [13] can all take place without in- 
volvement of metastable polymorphs. Further, simula- 
tions of charged colloids [14 show that sluggish dynam- 
ics can invalidate the closely-related Stranski-Totomanow 
conjecture [T5|, the prediction that the first phase seen is 
the one separated from the parent phase by the smallest 
free energy barrier. 

The limitations of these rules of thumb motivate us 
to look for connections between crystallization pathways 
and microscopic features of particle interactions and dy- 
namics. Here we study a lattice model of anisotropic 
particles. Our model is designed to mimic, in a generic 
way, the ability of materials such as proteins O |16] and 
ions [17] to crystallize by first forming a disordered phdise. 
It is also designed to be simple enough to allow thorough 
assessment of how its microscopic parameters control its 
crystallization behavior. Here we describe this behavior 
in detail. We show where in phase space Ostwald's rule 
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is likely to hold, and where it is likely to fail. Because 
the essence of the microscopic features of our model are 
common to a wide range of physical systems, the trends 
we have identifed might be used as a guide to select par- 
ticular crystallization pathways in simulation and exper- 
iment. 



II. MODEL AND SIMULATION PROTOCOL 

We consider a collection of particles that live on a fea- 
tureless two-dimensional substrate, which we model as 
a square lattice. Lattice sites may be vacant or be oc- 
cupied by a particle. Nearest-neighbor particles receive 
a 'nonspecific' interaction energy reward of —J. Parti- 
cles are anisotropic, and can point in any of R discrete 
directions. Neighboring particles receive an additional 
'specific' energy reward of —Q if they are aligned, and a 
penalty of -\-Q if they are antialigned. The larger is i?, 
the more precisely must two particles align before they 
receive the specific binding reward (we shall explore the 
effect of varying R). Each particle on the substrate feels 
a chemical potential — /i. 

We simulated our model using the Monte Carlo pro- 
cedure described in Appendix A. This procedure allows 
particles to translate (adsorb to and desorb from the sub- 
strate), and to rotate in place on the substrate. We have 
explored the effect of varying extensively the relative rate 
r of proposing rotation and translation moves, because 
we expect rotational and translational mobilities to vary 
considerably from one material to another. For instance, 
the limit of slow rotations (r <C 1) might be appropriate 
for particles, like DNA-linked colloids, that must unbind 
in order to rotate appreciably fi8\. It is also likely that 
particles of different sizes explore their positional and an- 
gular interaction ranges at different rates [T9 . We used 
this procedure in concert with umbrella sampling [20^ to 
calculate free energy landscapes for crystallization (bias- 
ing the size and degree of crystallinity of the system's 
largest cluster [8j, and in concert with transition path 
sampling [211 122] and forward flux sampling |23^ to gen- 
erate dynamical crystallization trajectories. The com- 
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bination of these methods reveals the distinct effect on 
crystahization of thermodynamics and dynamics. 



III. SIMULATION RESULTS 

In Fig. [l] we show the phase diagram of the model as 
a function of interaction strengths J and Q, for angular 
specificity R = 24. Snapshots of phases are also shown: 
particles with 4 parallel neighbors (crystalline particles) 
are green; particles with no parallel neighbors (fiuid par- 
ticles) are dark blue. Intermediate particles, with 1-3 
parallel neighbors, are blue-green. When both interac- 
tions are weak the stable phase is a homogeneous fiuid 
phase (H) of moderate density and little orientational or- 
der. When Q is large enough (above the 'freezing line'), 
the orient ationally-ordered solid (S) is stable. When J 
is sufficiently large (to the right of the 'demixing line'), 
phase H disappears and dense liquid (L) and sparse va- 
por (V) phases become viable. Our choice of fi (see Ap- 
pendix) ensures that above the freezing line and to the 
right of the demixing line the liquid phase lies intermedi- 
ate in free energy between the vapor phase and the stable 
solid. In this regime, Ostwald's rule suggests that a va- 
por should transform into a liquid before it crystallizes. 
But is this true? 

Crystallization pathways at points A, B and C in phase 
space are shown in Fig. [2] At point A in phase space, only 
phases H and S are viable, and crystallization is straight- 
forward: it consists of the direct transformation of H into 
S. At point B on the phase diagram (in our regime of in- 
terest) an empty substrate immediately becomes host to 
the metastable low-density vapor. We show in Fig. [2] the 
free energy of formation from the vapor of a dense nu- 
cleus as a function of nucleus size N and the number of 
crystalline particles Nc it contains. The minimum energy 
pathway ('mep', solid line) from the vapor to the crystal 
is an indirect one (sometimes called a 'two-step' or 'non- 
classical' pathway) that displays a liquidlike critical nu- 
cleus. The direct pathway (dashed line) via a crystalline 
critical nucleus is disfavored by about 15 /cbT. The indi- 
rect pathway is made possible by the intermediate liquid 
phase, but because we are far from the demixing line is 
not a result of critical density fiuctuations [24] . Dynam- 
ical trajectories generated using a wide range of particle 
rotation rates adhere to the indirect pathway (red line, 
fast rotation: r = 99; blue line, slow rotation: r = 0.01): 
a liquid nucleates on the substrate, and only subsequently 
does the crystal emerge from the liquid. For rapid rates 
of rotation the postcritical liquid readily transforms into 
a crystal while still only of small size [Ml] while for slug- 
gish rotation rates the liquid consumes the substrate, and 
fails to crystallize during the course of the simulation M2 , 

At point C in phase space the liquid remains interme- 
diate in free energy between the parent phase and the 
stable solid, but the driving force for crystallization is 
qualitatively different: the direct pathway, with a crys- 
talline critical nucleus, is preferred! The indirect path- 
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FIG. 1: Model phase diagram in the space of nonspecific (J) 
and specific (Q) attraction strengths, with phase snapshots, 
for fixed chemical potential (see Appendix). When both inter- 
actions are weak the stable phase is a homogeneous fiuid H of 
moderate density and no orientational order. 'Demixing' into 
a sparse vapor V and a dense liquid L (both orientationally 
disordered) can be induced by application of a strong non- 
specific attraction, J (the critical point on the horizontal axis 
is the regular Ising model one). The orientationally-ordered 
solid phase S is stable when the specific interaction Q is large 
enough, above the freezing line. We are interested in the re- 
gion of solid stability to the right of the demixing line. There, 
the sparse vapor phase lies above the dense fiuid phase in free 
energy, which in turn lies above the stable solid phase. Ost- 
wald's rule suggests that if we start from the vapor phase then 
the liquid will emerge prior to crystallization. But as we de- 
scribe below, this is not always true. The thermodynamically- 
favored critical nucleus is liquidlike below the horizontal red 
line, but crystalline above it. Points A, B, and C correspond 
to the simulations shown in Fig. [2] The connection between 
this phase diagram and the conventional temperature-density 
one is given at mean-field level in Fig. |A8| 



way with a liquidlike critical nucleus is still viable, but 
is disfavored by about bk^T. Because of this relatively 
small discrepancy in barrier heights, both pathways can 
be seen in dynamic simulations. For a sufficiently fast ro- 
tation rate (red line, r = 99) the direct pathway is taken: 
a crystal nucleates and grows on the substrate. No liquid 
is seen. For a slow rotation rate (blue line, r = 0.01) 
the indirect pathway is seen, and the substrate is again 
consumed by a liquid, [M3l |M4[ and | A7 



A relic of the liquid phase therefore inffuences the crys- 
tallization pathway some way past the freezing line, but 
eventually thermodynamics favors a direct mode of crys- 
tallization: the critical nucleus is crystalline above the 
horizontal red line on the phase diagram of Fig. [l] The 
question of how much liquid is seen in simulations is one 
that cannot be addressed by Ostwald's rule ('OR' in cap- 
tion). Assuming that it applies (it pertains to metastable 
intermediate phases, and the liquid at B and C is at most 
weakly so), it is satisfied in a macroscopic sense by the 
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FIG. 2: Crystallization pathways at points A, B and C in phase space of Fig. [T] At point A, only phases H and S are viable, 
and crystallization is straightforward: it consists of the direct transformation of H into S. Points B and C lie in the regime of 
interest, where the vapor, liquid and crystal phases lie in descending order of free energy. We show free energy surfaces (in a 
space of cluster size N and degree of crystallinity A^c) and dynamical trajectories for crystallization at phase points B and C, 
between which the thermodynamic mechanism for crystallization switches from an 'indirect' one to a 'direct' one. Snapshots 
bottom right show configurations generated using different dynamical protocols. The resulting trajectories sometimes uphold 
and sometimes violate Ostwald's rule (OR) and the Stranski-Totomanow (ST) conjecture. 



'slow rotation' (blue) trajectories at point B, but only in a 
microscopic sense by the 'fast rotation' (red) trajectories. 
It is violated by the 'fast' trajectories at C, but not by the 
'slow' trajectories. Moreover, the latter involve passage 
over a free energy barrier larger than the smallest avail- 
able, going against the sense of the Stranski-Totomanow 
conjecture ('ST' in caption). 

The microscopic thermodynamic control of our model's 
crystallization pathway is the competition between the 
angular specificity R and potency Q of specific binding. 
In Fig. [3] we show that the thermodynamically-preferred 
composition of the critical nucleus (and of small nuclei 
of specified size) is liquidlike for small Q (just above 
the freezing line), and becomes increasingly crystalline 
as one increases Q: this trend acts against the sense of 



Ostwald's rule. By contrast, increasing R acts in favor 
of Ostwald's rule: as R grows, the liquid becomes more 
strongly metastable to crystallization, giving rise to free 
energy barriers that hinder the transformation of a liquid 
droplet into a solid one (compare surfaces top right and 
top left). 

IV. ANALYTIC RESULTS 



Interestingly, we can anticipate these thermodynamic 
trends using simple microscopic theory. We calculated 
analytically (see Appendix B) the model's bulk free en- 
ergy /eff, in a self-consistent mean- field approximation, 
as a function of order parameters p (density) and r (crys- 
tallinity): 



fes{p,T) = ^ (Jzp^ + QzT^) - fceTln (1 + e^^''P+^^^ [i? - 2 + 2cosh(/?Q^r)]) . (1) 



Here 2: = 4 is the lattice coordination number. From which resembles qualitatively its simulated counterpart. 
Eq. ([1]) we obtained the phase diagram shown in Fig. |4j 
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FIG. 3: Microscopic thermodynamic controls of crystallization. Thermodynamically-preferred composition (/b, left, is fraction 
of crystalline interactions, and A^c, right, is number of crystalline particles) of a critical droplet A^crit or one of specified size 
A/", as a function of specific interaction strength Q, for three values of angular specificity R. The bulk solid is stable to the 
right of each freezing line (vertical gray dots labeled 'FL'), but a droplet of finite size prefers to be liquid some way past the 
phase boundary. As Q increases, however, small droplets prefer to be solid, a trend that opposes Ostwald's rule. By contrast, 
increasing R acts in favor of Ostwald's rule by rendering the liquid more strongly metastable to crystallization: corresponding 
free energy surfaces (at right) possess barriers that hinder the crystallization of a liquid droplet. 



The microscopic parameters of the model control dis- 
tinct critical behaviors, but they also shape the bulk free 
energy landscape even in regimes away from any phase 
transition. Furthermore, previous work reveals that the 
bulk free energy landscape suggests the qualitative char- 
acter of the thermodynamically-preferred crystallization 
pathway [25l [26]: bulk wells indicate driving forces for 
appearance of phases, and bulk barriers are an impor- 
tant component of droplet surface tensions in a Cahn- 
Hilliard approximation. On bulk surfaces we calculated 
the minimum energy pathways between parent and solid 
phases (shown by white dashed lines on surfaces at points 
A, B, C), and shaded the phase diagram according to 
how 'direct' are those pathways. The resulting 'path- 
way diagram' describes qualitatively the thermodynamic 
driving force for crystallization for given microscopic pa- 
rameters, and mirrors the changes of cluster composition 
seen in our umbrella sampling (thermodynamic) simu- 
lations: indirect crystallization pathways become viable 
to the right of the demixing line, and are supplanted 
by a direct mechanism some distance above the freezing 
line. This correspondence is a central result of this pa- 
per. This change of mechanism takes place in a regime 
of phase space in which the hierarchy of stable phase re- 
mains unchanged: in other words, the phase diagram in 
this regime is featureless, but the pathway diagram is not. 

The analytic theory allows a comprehensive survey of 
parameter space. Increasing i?, the number of accessi- 
ble orientational states, has the effect of increasing the 
bulk free energy barrier between crystal and liquid phases 
(Fig.lsl, a trend reflected in simulations of finite-size clus- 
ters (Fig. [3|. Fig. [g] shows pathway diagrams for three 
particles with different value of which display similar 
trends (though a stronger specific interaction is neede d to 
effect a change of pathway for large R). Finally, Fig. A8 
shows the relation between the J, Q phase diagram and 



the conventional temperature-density one. The analytic 
estimate of pathway is a useful starting point for directing 
simulations, but it contains no dynamical information: 
we have seen in simulations that sluggish rotational dy- 
namics can render the actual crystallization mechanism 
different to the thermodynamically-preferred one. 



V. CONCLUSIONS 

We have identified the microscopic controls of crys- 
tallization in a model of anisotropic particles that can 
form both disordered and crystalline phases. Our work 
complements previous studies that reveal how changes of 
intensive parameters, such as temperature [27] and pres- 
sure [28 , can change crystallization pathways in model 
systems. Although rules of thumb are of limited use in 
predicting our model's crystallization pathway, a combi- 
nation of simulation and analytic theory reveals that a 
liquid phase is likely to be seen prior to crystallization 
if 1) particles rotate sluggishly; 2) if particles must align 
precisely in order to crystallize (i.e. if R is large); and 
3) if the specific attraction Q is just strong enough to 
render the crystal stable (but not so strong as to render 
small droplets of crystal lower in free energy than droplets 
of liquid of similar size). While our model is idealized, 
and certainly does not capture the detailed microscop- 
ies of the way real particles interact, it does show that 
the type of complex crystallization pathway observed in 
many experiments emerges as soon as one assigns to a 
particle a translational degree of freedom and a rota- 
tional one. Further, while we expect differences in the 
details of particle interactions in two- and three dimen- 
sions, the correspondence of mean field theory (which 
tends to work better the higher the dimension) and sim- 
ulation results, and related dynamic 'crossovers' seen in 



distance 



distance 



distance 



□ direct 



indirect 










r 







0.2 0.4 o.e 

p 




0.2 0.4 0.6 0.8 1 



FIG. 4: Visualizing the shaping of crystaUization landscapes by microscopic parameters. Mean-field 'pathway diagram' (bottom 
left, R — 24) and free energy surfaces at points A, B, and C in phase space (mep is dotted white line; star denotes transition 
state), calculated using Eq. ([T]). Above each surface we show free energy along the mep (solid black line with star), and, for 
comparison, free energy along a direct pathway (B) and an indirect pathway (C) (shown dotted; these pathways are not shown 
in panels below). Pathway diagram shading indicates the preferred crystallization mechanism; the plot above it shows where 
on the lines A ^ B ^ C the transition state rtr is crystalline. The qualitative trends identified mirror those seen in our 
simulations. 



three-dimensional models [27], suggest that the results 
of this study are not limited to two-dimensional systems. 
With this in mind, we conjecture that the observations 
made here might be used as a starting point to guide 
experiments. For instance, they suggest that one could 
control the crystallization pathway of a protein if one 
could devise ways of altering, independently, that pro- 
tein's nonspecific attraction 'J' (e.g. using PEG) and 
specific attraction 'Q' and 'i?' (e.g. through mutation 
or using multivalent salt 129 ). One might also design 
anisotropic nanoparticles of specific rotational and trans- 
lational mobilities in order to encourage a direct or indi- 
rect crystallization pathway. Finally, we used the 'path- 
way diagram' of Fig. |4] to select interesting regions of 
phase space for our simulation work: analytic treatment 
of other models might furnish similar microscopic 'maps' 
to guide simulation studies. 
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FIG. 5: Effect of varying i?, the number of accessible rota- 
tional states. We plot free energy (vertical) along the min- 
imum energy pathway ('distance' axis) of a mean-field free 
energy surface at J = 3kBT, for different values of R. The 
value of Q is such that we sit just above the freezing line for 
each R (see Fig. [I] for = 24). In all cases the pathway 
observed is an indirect one from vapor (V) to liquid (L) to 
solid (S). The first barrier seen corresponds to the vapor-to- 
liquid transformation, and the second to the liquid-to-crystal 
transformation. Increasing R leads to the growth of a large 
entropic barrier to crystallization. 
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J/UbT J/UbT J/UbT Q/UbT 

FIG. 6: Summarizing the shaping of crystaUization landscapes by microscopic parameters. We show mean- field phase diagrams 
for three values of i?, shaded according to where direct (light) and indirect pathways (dark) are favored. Right panel: barrier 
height AC between liquid and crystal phases, as a function of Q, for J = Ak^T (vertical lines on the three diagrams). The 
change of crystallization driving force from direct to indirect can both preempt (large R) and follow (small R) the loss of liquid 
metastability. 
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Appendix A: Simulation details 

We simulated our model using the following grand 
canonical Metropolis Monte Carlo procedure. This pro- 
cedure effects a diffusive dynamics and assumes the sub- 
strate to be in contact with a thermal bath and a par- 
ticle bath (i.e we assume the substrate to be in con- 
tact with bulk solution). We select at random a lat- 
tice site. If that site is occupied by a particle then 
with probability pdei < 1 we attempt to delete its occu- 
pant. We accept this deletion with probability Pdeiete = 
min(l,^exp(-/3A£;-/3/i)), where /3 = l/{k^T) 
and is the change of interaction energy following 
the proposed deletion. With probability 1 — pdei we 
instead attempt to change the particle's orientation by 
±1 unit, modulo R (we assume particles to rotate in 
a plane, and so orientation R neighbors orientation 1). 
We accept changes of rotation with probability Protate = 
min (1, exp {—PAE)). If the chosen site is vacant then we 
attempt to occupy it with a particle whose orientation is 
chosen randomly. This attempt succeeds with probability 
^insert = miu (1, pdei^ Gxp (— /3A£^ + . The factors 
of pdei ill insertion and deletion rates are required to pre- 
serve detailed balance: insertions are always attempted if 
a lattice site is vacant, but deletions are attempted only 
with probability pdei < 1 if a lattice site is occupied. The 
factors of R are present for a similar reason: insertion of 
a particle of a given orientation is attempted with prob- 
ability 1/R^ but proposing the reverse of that particular 
insertion occurs with unit probability. We used a lat- 
tice of = (100)^ sites, periodically replicated in each 
direction. 

The basic rate for particle translations (adsorptions 
and desorptions) is ~ pdei- The basic rate for parti- 
cle rotations has two contributions: the first scales as 
(1 —pde\)/R^i and comes from explicit rotation moves; 
the factor of R^ accounts for the characteristic time to 
visit R rotational states. The second contribution scales 
as Pde\/R^ and comes from explicit translations (parti- 
cles attach to the substrate with randomly-chosen orien- 
tations, and so removal and reattachment of a particle 
allows an effective sampling of its orientation). The first 
mode of orientation-sampling is most effective in the bulk 
of a liquidlike cluster, while the second mode operates 
most readily at the surface of a cluster (where particle 
detachments are most frequent). There is therefore no 
constant effective rate at which a cluster explores the 
configuration space (TV, TVc) shown in Fig. 1. We report 
values of a parameter r = (1 — PdeO/Pdeh the relative 
rate of proposing a rotation or translation. Our scheme 
does not account for in-plane particle diffusion, focusing 
instead on mass transport from the bulk. 

We set /i = —Jz/2 — k^ThiR^ where z = 4 is the coor- 
dination number of the lattice. Our choice of /i ensures 
that, to a mean- field approximation, liquid and vapor 
phases are equal in free energy. The contribution —Jz/2 
is the usual Ising model term. The term —k^ThiR pe- 



nalizes particles relative to vacancies by an amount that 
exactly compensates the entropy difference between par- 
ticles and vacancies. This choice is motivated by the 
fact that we consider our simulation protocol to reflect 
the diffusion of material to and from the substrate, with 
no change of that material's rotational freedom, rather 
than to model its creation or destruction (which would 
be accompanied by creation or destruction of rotational 
entropy). In simulations, orientational correlations in the 
liquid lower its free energy below that of the vapor by an 
amount that decreases with increasing R. 

In umbrella sampling simulations, cluster size N and 
crystallinity A^c were constrained using harmonic bonds 
of spring constant 0.2 within windows of width 5 (par- 
ticles). Dynamical trajectories were sown, grown and 
harvested using the original forward flux sampling (FFS) 
algorithm [30 and the aimless shooting transition path 
sampling algorithm [311 [32]. FFS simulations used an 
interval of 5 between interfaces, 1000 starting configura- 
tions in the initial basin and 10^ trials per interface. 

Appendix B: Mean-field theory 

The energy function of our model is H = 
Zlili {\ Y^j Uij - IJin^, where j runs over the z = 4 
nearest neighbors of i, and /i is a chemical potential. The 
variable is if lattice site i is vacant, and is 1 if it is 
occupied. The pairwise interaction Uij is 

Uij = -niUj (^J + QS{si, Sj) - Q5{si, sj)^ , 

where Si = 1, 2, is the orientation of the particle at 
lattice site i. The function S{si^ Sj) is 1 if 5^ = Sj (aligned 
particles receive an extra energetic reward), and is zero 
otherwise. The function 6{si, Sj) is 1 if Si and Sj are R/2 
units different, modulo R (antialigned particles receive 
an energetic penalty), and is zero otherwise. 

We can derive the free energy of this model in a mean- 
field approximation [33l IM] by assuming that each site 
feels only the thermal average of the fluctuating variables 
at neighboring sites. The effective interaction at a given 
site is to this approximation 

/7eff = — Jzn{n) 

R 

- Qz^n5{s,q) (^{n5{s,q)) - {n5{s,q))^ , 

where we have dropped site labels. By symmetry, all 
but two terms in the sum over q vanish. The ef- 
fective free energy per site is then /eff = E — TS^ 
where E = ^(/7eff) — /J^{n) and —TS = /cBT(lnPeq), 
where Peq = e'^'^-'' /Tt e'^^-^' and Heff = t/eff - /in. 
Thermal averages are defined self-consistently through 
the relation (A) = Tr(APeq). The trace Tr(-) = 

5^n=o,i Z)f=i +^n,o} (•) can be carried out by as- 
suming, without loss of generality, the ordering direction 
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FIG. A7: Distributions of the number of crystalline particles in a cluster, A^c, from dynamical trajectories at phase point C 
(see Fig. 1, main text). Distributions correspond to clusters at criticality, A^crit, and twice the critical size, 2A^crit. For sluggish 
rotation rates, distributions are peaked at low values of Nc'. on average the growing nucleus is liquidlike. For sufficiently rapid 
rotation rate the direct pathway is preferred. In all cases, trajectories can be observed that buck the trends shown (note the 
tails of each distribution). 
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FIG. A8: Supplement to Fig. |4] relation of the J, Q phase diagram to the more conventional T, p one. For two different 
varying-temperature lines {J/Q = 2.3 and 0.8) along our mean-field phase diagram (left, R — 24), we show phase diagrams 
in the density (p)-temperature (T) plane. These phase diagrams resemble qualitatively those of a protein (middle) or argon 
(right), with their respective stable- and metastable fluid-fluid demixing critical points [35]. We see that depending on where 
one lies in parameter space, changing temperature can change the thermodynamically preferred crystallization pathway in a 
complicated way. 



to be 5 = 1. The result is Eq. ([T]), in which p = (n) is Appendix C: Supplementary Figures 

the density and r = {nS{s^ 1)) — {nS{s^ 1)) is the crys- 
tallinity order parameter. The latter distinguishes dis- 
ordered fluid phases (for which r = 0) from the ordered 
solid phase (for which r 7^ 0). 



